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Abstract 

We demonstrate the nontrivial scaling behavior of Ising models defined on (i) a 
donut-shaped surface and (ii) a curved surface with a constant negative curvature. 
By performing Monte Carlo simulations, we find that the former model has two 
distinct critical temperatures at which both the specific heat C(T) and magnetic 
susceptibility x(T) show sharp peaks. The critical exponents associated with the 
two critical temperatures are evaluated by the finite-size scaling analysis; the result 
reveals that the values of these exponents vary depending on the temperature range 
under consideration. In the case of the latter model, it is found that static and 
dynamic critical exponents deviate from those of the Ising model on a flat plane; 
this is a direct consequence of the constant negative curvature of the underlying 
surface. 
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1 Introduction 



The recent progress in fine processing technologies has enabled the fabrica- 
tion of nanoscale materials with novel shapes [1,2,3,4,5,6,7,8,9]. A study of the 
basic properties of these materials is crucially important for many potential 
applications; it also helps to gain a new perspective on nanodevice modeling. 
From a fundamental point of view, it is interesting to determine how the geom- 
etry of these nanostructures influences their physical properties, particularly 
their critical behavior close to the phase transition. It is known that at the 
phase transition, certain physical quantities exhibit power-law behavior that is 
universal depending on certain essential symmetries of the system [10,11,12]. 
Hence, an alteration in the geometric symmetry of the system possibly results 
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in an alteration in critical behavior. On this background, several studies have 
suggested that systems with novel geometry exhibit peculiar critical behav- 
ior [13,14,15,16,17,18,19,20,21]. Nevertheless, no systematic investigation has 
been conducted on the relationship between geometry and critical phenomena 
that occur in magnetic nanostructures. 

In this context, one of the most interesting issue is to study the critical prop- 
erties of low-dimensional magnetic systems with nonzero geometric curvature. 
The two-dimensional Ising model with ferromagnetic interaction is a typical 
example of a system capable of a magnetic phase transition. It is well known 
that this model undergoes the ferromagnetic transition in zero external field 
as the temperature decreases below a critical temperature T c [22,23]. In the 
vicinity of T c , thermodynamic quantities exhibit the power-law behavior that 
is characterized by critical exponents. The values of these exponents are deter- 
mined only by certain global symmetries (and dimensionality) of the system; 
this has been confirmed both theoretically and experimentally in the case of 
the Ising model on a flat plane [24,25,26,27]. Currently, we are interested in 
the effect of geometry on the critical exponents, that is, whether a geomet- 
ric alteration induces a quantitative alteration in the values of the critical 
exponents of the system. 

In this study, we investigate the critical behavior of the Ising models defined 
on two types of curved surfaces: (i) a donut-shaped surface produced by glu- 
ing a flat plane and (ii) a surface with a constant negative curvature (i.e., 
a pseudosphere). Monte Carlo simulations combined with finite-size scaling 
analysis reveal that in both Ising models, the value of the critical exponent 
7 of the zero-field magnetic susceptibility x(T) (as well as the dynamic criti- 
cal exponent z) deviates from that of the conventional Ising model on a flat 
plane. Furthermore, the donut-shaped Ising lattice model has two distinct crit- 
ical temperatures at which the specific heat and magnetic susceptibility show 
sharp peaks. These surprising results indicate that the geometry of the under- 
lying surface is significantly related to the phase transition of low-dimensional 
magnetic materials. 



2 Donut-shaped Ising models 

First, we consider a square Ising lattice model defined on a donut-shaped 
surface (see Fig. 1) [28]. The donut-shaped lattice is constructed by gluing 
the edges of a square lattice on a flat plane; hence, this lattice is topologi- 
cally equivalent to a square flat lattice with the periodic boundary condition. 
However, this equivalence breaks down when the difference in the intrinsic ge- 
ometry between the flat plane and the donut-shaped surface is considered. The 
gluing procedure does not preserve the separation of the neighboring sites 
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Fig. 1. Schematic illustration of the square lattice defined on a donut-shaped sur- 
face obtained by gluing the edges of the flat plane. The spacial modulation of the 
separation occurs only along the X direction. 



i and j, thereby resulting in the spatial modulation of the coupling strength 
Jjj of the nearest-neighbor spins Sj and Sj. To our knowledge, it has not yet 
been clarified how the gradual modulation of affects the thermodynamic 
properties of the embedded Ising model. (In contrast, the Ising models con- 
sisting of random or alternative bonding have been extensively studied thus 



where the constants a and Jo represent the separation of neighboring spins 
and the coupling constant of the square lattice before the gluing procedure, 
respectively (in other words, these constants belong to the square lattice on 
the flat plane). In the definition of J^, the effect of the spatial modulation 
of Tij is taken into account so that Jy becomes spatially dependent. In the 
actual calculation, we have constructed donut-shaped Ising lattices with dif- 
ferent sizes by gluing planar square lattices with linear dimensions of L = 200, 
250, and 300 in units of a. For all donut-shaped lattices, the minimum and 
maximum separations are fixed as r min = 0.88 and r max = 2.88, respectively 
(see Fig. 1). 



far [29,30,31].) 



The donut-shaped Ising model is described by the Hamiltonian 



with Jaa = 

r ■ ■ 
1 i J 
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Fig. 2. (a) Specific heats and (b) magnetic susceptibilities of the donut Ising model. 
There appear two sharp peaks at T c \ = 1.47 and T C 2 = 2.40. 

3 Occurrence of two critical temperatures 



We have performed Monte Carlo simulations [32] to evaluate the specific heats 
C{T) and magnetic susceptibilities x(T) defined by 



C(T) 



(E 2 ) - (El 



(2) 



and 



X(T) 



(m 2 ) - (My 



(3) 



Here E and M are the internal energy and the spontaneous magnetization of 
the system, respectively; the angular brackets indicate to take thermal average. 
Figure 2 shows the numerical results for various system size; L=200 (open 
triangles), 250 (solid circles), and 300 (open squares). Most striking is the fact 
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Fig. 3. Spatial configurations of Ising spins on the donut surface. Black and white 
dots represent the positive and negative sign of the Ising spin, respectively. Each axis 
corresponds to the direction of X and Y on the donut surface depicted in Fig. 1. 
The horizontal lines Y = and Y = ±150 represent the internal and external 
circumference of the donut surface. 

that both data exhibit two sharp peaks at T c \ = 1.47 and T C 2 = 2.40 in unit 
of Jo/kB- The values of these critical temperatures are well apart from that 
of the square Ising model on a flat plane [32] : 



T, 



ln(l + y/2) 



2.27. 



(4) 



We have also attempted to evaluate the critical exponents of x{T) by means 
of the finite-size scaling analysis. The scaling hypothesis states that in the 
vicinity of T c , x for a finite system size L obeys the scaling form [10]: 

X (T,L)^L^f(\T-T c \L^) } (5) 



where 7 and v are referred to as the critical exponent of xCO an d that of the 
correlation length £(T) of the order parameter, respectively. Intriguingly, in 
our donut Ising model, the value of the exponent 7 depends strongly on the 
temperature range. For instance, at immediately lower and higher than T&, 
the exponent yields 7 ~ 3.0 and 7 ~ 0.6, respectively. On the other hand, 
it yields 7 ~ 0.9 at just higher than T c2 . These values of 7 obviously deviate 
from that for the planar Ising model: 72^ = 7/4. 

The occurrence of the two distinct critical temperatures in the donut-shaped 
Ising model is qualitatively explained by observing the spatial configurations 
of the Ising spins S{ at different temperatures. Figure 3 shows the spatial 
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configurations of the Ising spins on the donut surface; black and white dots 
represent the positive and negative signs of the Ising spin, respectively. At 
T > T c2 (in the right panel), the spin directions are completely random and 
therefore the system has no spontaneous magnetization. On the contrary, at 
T < T c i (in the left panel), most spins are oriented in the same direction; 
hence, the system has long-range order. The configuration in the intermediate 
temperature range of T c \ < T < T C 2 (in the middle panel) is important. In the 
last large ferromagnetic domain exists along the internal circumference 

of the donut-shaped surface (Y ~ 0), while the region close to the external 
circumference (Y ~ 150 or —150) remains in a paramagnetic state. 

The partial growth of the ferromagnetic domain results from the gradual mod- 
ulation of Jij over the donut-shaped surface. Close to the internal circumfer- 
ence, the magnitude of is sufficiently large to overcome the thermal fluc- 
tuation; hence, the ferromagnetic domain at the inner part persists until the 
temperature reaches T c2 . On the other hand, the ferromagnetic order along 
the external circumference disappears at T c i, since in the latter region is 
very weak. Consequently, the two critical temperatures are produced by local 
ferromagnetic transitions near the internal and external circumferences of the 
donut-shaped surface. 

It is emphasized that the peculiar behaviors of the donut-shaped Ising model, 
such as the two peaks of x(T) and the quantitative shift of 7, are direct con- 
sequences of the gluing procedure; in other words, they are caused by the 
mechanical deformation of a flat magnetic plane. Hence, these findings imply 
that it is possible to control the physical properties of magnetic nanostruc- 
tures by inducing a local (or global) mechanical deformation. The relevance 
of our results to actual magnetic systems with novel geometry is also under 
consideration. 



4 Ising models on negatively curved surfaces 

In the previous sections, we have considered the donut-shaped Ising lattice 
model in which the separation of neighboring spins is spatially dependent 
due to the external gluing procedure. Alternatively, however, it is possible to 
construct the regular Ising lattice model on certain curved surfaces, wherein 
the Ising spins are assigned with an equiseparation over the entire surface. A 
typical example is an Ising model defined on negatively curved surfaces, i.e., 
on a pseudosphere [33,34]. The pseudosphere is a simply connected infinite 
surface whose Gaussian curvature at any point has a constant negative value 
[35,36]. More importantly, in the latter model, the surface curvature effects 
manifest themselves in scaling behavior in a different way from that in the 
donut-shaped Ising model. In what follows, we will give a brief review of our 
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Fig. 4. Local bonding structures of regular (a) heptagonal and (b) triangular lattices 
in terms of the Poincare disk representation. The size of our lattice is determined 
by the number of concentric layers of polygons, denoted by r. 

previous study on the Ising model on the pseudosphere [37,38]. By comparing 
the results obtained for the two distinct curved surfaces, we gain a deeper 
understanding of the effect of geometry on the Ising model. 

Figure 4 illustrates (a) regular heptagonal and (b) regular triangular Ising 
lattices embedded on the pseudosphere in terms of the Poincare disk represen- 
tation [39]. Although polygons depicted in the figures appear to be distorted, 
they all are surely congruent in the sense of the intrinsic geometry on the 
pseudosphere. Onto these lattices, we assign the ferromagnetic Ising model 
with nearest neighbor interaction, described by the Hamiltonian 



and perform Monte Carlo simulations [40] to evaluate the critical exponent 
7 of the magnetic susceptibility x- It is noted that in the present case, the 
scaling form of x should be expressed as 



where \i describe the divergence of the correlation volume as £v(T) oc |T — 
T C |~ M . Similar to the case of the donut Ising model, quantitative evaluation 
of these exponents can be achieved by using the finite-size scaling technique 



Before addressing numerical results, it deserves comment on the boundary 
contribution in the Ising model on the pseudosphere. It follows from Fig. 4 
that the size of our lattice is determined by the number of concentric layers of 
polygons, denoted by r. When r ^> 1, the number of total sites is approximated 
as N(r) oc e r ; this exponential increase in N means that the contribution from 
the boundaries survives even in the thermodynamic limit r = oo. In order to 



(6) 




(7) 



[37,38,40]. 
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Fig. 5. <5r-dependences of (a) 7 and fi and (b) z. Each data point was extracted by 
means of the finite size scaling analysis for the system sizes 4 < rj n < 8. 

extract the bulk critical properties, therefore, we have taken into account only 
the Ising spins involved in the interior r in layers (r in < r out = r in + 5r) when 
performing the scaling analysis. Asymptotic behaviors of 7 and [i for large Sr 
(as well as large r- m ) provide estimations of the bulk critical exponents. 

Figure 5 (a) shows the values of critical exponents [i and 7 as a function of the 
number of disregarded layers Sr. We see that, while the curve of /1 converge 
to a particular value of fi ~ 2 or less, that of 7 has no tendency to converge 
to 72d = 7/4 for large 5r. Notably, these results are consistent in quality 
with the conjecture deduced by the quantum-field theory [41]. It states that 
the Ising model embedded on the pseudosphere should yield the mean-field 
critical exponents (i.e., 7mf = 1 and /imf = 2) when the boundary contribution 
may be omitted. These mean-field behaviors can be simply attributed to the 
relation N cc N s between N and the number of spins along the boundary N s 
0, which implies that the pseudosphere act as an infinite-dimensional surface. 

It is also interesting to evaluate the dynamic critical exponent z of the Ising 
model on the pseudosphere, which can be done by using the short-time re- 
laxation method [42,43]. The 8r- dependence of z are summarized in Fig. 5 
(b). We observe that z monotonically decreases with 5r; this indicates that 
for sufficiently large 5r, z takes a value considerably smaller than that of the 
planar Ising models: z = z/d ~ 1.1. That is, z on the pseudosphere is also 
expected to exhibit the mean field exponent zmf = = 1/2 at 5r 1, 

since zmf = 2 for the Ising model and d c = 4 the upper critical dimension. 
Quantitative determination of z for 5r ^> 1 as well as those of other static 
exponents (such as j3 and rj) are under consideration; preliminary results are 



1 Since the relation N s oc iV 1 ^' d ) holds for an Ising lattice in d dimension, N s oc N 
effectively consequences d = 00. 
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given in Ref. [44]. 



5 Conclusion 

To summarize, we have investigated the critical behaviors of the two-dimensional 
Ising models embedded on a donut-shaped surface and a pseudosphere. In the 
case of the donut-shaped Ising model, the specific heat C(T) and magnetic 
susceptibility x(T) exhibit two sharp peaks at T cl = 1.47 and T c2 = 2.40. 
The two critical temperatures are produced by the spatial modulation of the 
coupling strength that leads to the formation of a partial ferromagnetic 
domain along the internal circumference of the donut-shaped surface. Further- 
more, we have attempted to evaluate the critical exponent 7 of x(T). When 
the temperature is less than T c i, the exponent yields 7 ~ 3.0; however, when 
the temperature is greater than T cl , it yields 7 ~ 0.6; further, at a tempera- 
ture just above T c2 , it yields 7 ~ 0.9. Evidently, these values of 7 deviate from 
those of the planar Ising model. We have also demonstrated that in the case of 
the Ising model on the pseudosphere, the static exponent 7 and the dynamic 
exponent z deviate from those of the planar Ising model. We believe that the 
results of our study provide new information for the theoretical development 
of critical phenomena that occur on general curved surfaces. 
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